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We propose to correlate transmittance maps and spectral-density maps of planar junctions, in 
order to analyze quantitatively and in detail spin-dependent transport calculations. Since spectral- 
density maps can be resolved with respect to atom, angular momentum, and spin, the resulting 
correlation coefficients reveal unequivocally, e. g., which layers or which orbitals determine the tunnel 
conductances. Our method can be used for transport calculations within the Landauer-Biittiker 
formalism. Its properties and features will be discussed by means of a pure bcc Fe(OOl) lead as well 
as an extensively studied Fe(001)/MgO/Fe(001) planar tunnel junction. 
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I. MOTIVATION 

Spin electronics — or spintronics for short — is one of the 
major topics in contemporary physics (see for example 
Refs. |T]and|2)). With respect to both device applications 
and fundamental physics, planar junctions have been and 
are being investigated experimentally and theoretically 
with great effort. The spin-dependent transport prop- 
erties of magnetic tunnel junctions show up as tunnel 
magnetoresistance (TMR), that is the change of the con- 
ductance upon reversal of the magnetization direction 
in one of the two electrodes^ ^. Replacing the insulat- 
ing barrier by a ferroelectric, a tunnel electroresistance 
(TER) effect can be observed in addition. In this case, 
the conductance depends as well on the orientation of the 
electrical polarization in the ferroelectric^. 

In TMR experiments, the current that is flowing 
through a tunnel device is detected in dependence on ex- 
ternal fields and device parameters (e.g. bias voltage and 
individual layer thicknesses). The current- voltage char- 
acteristics of tunnel devices are often interpreted within 
the Julliere modeP, perhaps due to its simplicity. The 
validity of this model has been severely questioned'' be- 
cause it relates the TMR ratio exclusively with the spin 
polarization of the electrodes and, thus, neglects the in- 
terface and barrier regions completely. If the Julliere 
model would be valid, the tunnel conductance would not 
depend on the interface material at all, in contrast to ob- 
servations. In other words, the interface region is essen- 
tial and has to be described in theory as good as possible. 

From the preceding it is evident that a reliable de- 
scription of transport in planar junctions must capture 
the essential properties on an atomistic level. Hence, 
present theoretical investigations of transport rely on so- 
phisticated first-principles approaches to the electronic 
and magnetic structures. 

Advanced first-principles approaches to transport al- 
low a very detailed analysis of the electronic structure 
and the conductance. In particular, the probably most 
important question — which orbitals in which layer deter- 



mine the transport properties — can be answered. How- 
ever, the amount of output data that is produced by 
modern computer codes becomes often unmanageable for 
human beings. Therefore, one restricts oneself to repre- 
sentative subsets or comprises the numerical data into 
manageable representations. For example, one puts the 
focus of the analysis on a single wavevector in the two- 
dimensional Brillouin zone (2BZ), typically to the 2BZ 
centei^. By doing so, one should be aware that such a re- 
striction could fail because considerable parts of the 2BZ 
may contribute to the conductance. Examples for data 
representations are transmittance and spectral-density 
maps in the 2BZ (introduced in Section [ll|) which have 
become established analysis tools. They allow in prin- 
ciple to answer the above question. Because they are 
compared visually, they leave space for speculation and 
interpretation; or phrased differently, they introduce am- 
biguity. Apparently, there is need for an improved analy- 
sis tool which allows to analyze transport properties un- 
equivocally and quantitatively, rather than ambiguously 
and qualitatively. 

In this paper, we propose a quantitative analysis of 
transport properties which goes beyond the approaches 
sketched above. We propose to correlate transmittance 
and spectral-density maps. The resulting correlation co- 
efficients reveal unequivocally, for example, which lay- 
ers or which orbitals determine the tunnel conductances. 
Our method can be used for any computer code which 
relies on a Landauer-Biittiker-type approach, thus being 
applicable in most of the present-day transport calcula- 
tions. Its properties will be discussed for two junctions 
exhibiting a planar geometry. 

The paper is organized as follows. In Section [TT] we 
introduce our approach for analyzing transport calcu- 
gives a brief overview of our nu- 



lations. Section HI 



merical approach. The correlation analysis is applied 
in Section [lV| to a pure Fe(OOl) lead (IV A) and to an 
Fe(001 )/MgO/Fe(001) magnetic tunnel junction (MTJ) 
(IV B). Concluding remarks are given in Section [V[ 
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II. THEORETICAL ANALYSIS OF THE 
TUNNEL CONDUCTANCE 

We consider a planar junction which consists of a left 
electrode £, an interface region X, and a right electrode 
IZ. Due to the translational invariance parallel to the 
interface region, the Bloch states in the electrodes are 
indexed by the (in-plane) wavevector k\\ in the 2BZ; k\\ 
is conserved in the scattering process. 

The bias voltage V opens an 'energy window of trans- 
port'; the chemical potentials jic and fin of >C and 7^, 
respectively, differ by eV = i^c — l-^n- Without loss of 
generality we consider the case V > 0, for which incom- 
ing occupied Bloch states in C can be transmitted into 
outgoing unoccupied Bloch states in IZ. 

According to Landauer and Biittikei^^, the conduc- 
tance C is given by 



C{V) 



{V]E,k\\)d^kdE. (1) 



J2BZ 



The wavevector integral is over the two-dimensional Bril- 
louin zone. The transmittance T^^^(V; /cy) is the 
sum over the transmission probabilities of all incoming 
occupied states A in £ and outgoing unoccupied states 
p in IZ. It is related to the scattering matrix S of the 
interface region by 

T^-^(F; k\\) = Yl I^aV(^; ^ll) r • (2) 

Xp 

The above transmittance is a key quantity in the pro- 
posed analysis; it is conveniently displayed versus k\\ at 
fixed E and V ^ in so-called transmittance maps (T maps). 
Note that T can be regarded as a global quantity since it 
depends on the electronic structure of the entire junction, 
that is both electrodes (incoming and outgoing Bloch 
states) and the interface region (scattering matrix). 

The local electronic structure of the device is described 
in terms of the spectral density 



Nc{V;EM) 



1 



ImTr«G+(y;E,fc||). 



(3) 



a is a compound index which can comprise for example 
layer, atom, orbital, angular momentum, spin indices or 
point-group representation. For a given division of the 
complete index space, the subsets of a are disjoint. The 
trace Trc^ of the Green function is restricted to the 
given a and, thus, allows a detailed analysis of the elec- 
tronic structure. As for the transmittance, the spectral 
density is displayed in maps versus k\\ {N maps). 

In previous publications we have analyzed the conduc- 
tance by relating pronounced features in a transmittance 
map with features in the associated spectral-density 
maps. If these features show up in both the T map and 
in an Na map we concluded that the chosen set a deter- 
mines the transmittance in this (V, E) region. This way 
we could show that for Fe/Mn/vacuum/Fe junctions the 



topmost Mn layer governs their TMR ratio'^. However, 
for more complicated systems it turned out that the vi- 
sual inspection of the maps becomes ambiguous, tedious, 
and not very reliable. As a consequence, we propose to 
correlate properly normalized T and N maps by means 
of a projection (inner product). This procedure results 
in a set of a few unambiguous numbers which allow to 
determine rapidly the set of significant a indices. 

For given energy E and bias voltage V we define the 
average value of a function X{V; E^k\\) of the transmit- 
tance T(V; E^k\\) and the spectral density Na{V; E,k\\) 
as the average over the 2BZ, 



^^2BZ J2BZ 



(4) 



where 1^2BZ is the area of the 2BZ^^. The correlation 
coefficient Ca{V;E) is then defined by 



Ca{V;E) 



(5) 



This quantity can be interpreted as an inner product of 
the normalized TiV] ^•>'^\\) and No^iV] E^k\\). Since the 
latter are semi-positive for all /cy, < Ca{V;E) < 1. 

Note that Ca{V; E) is invariant with respect to scaling 
T{V]E,k\\) and No,{V]E,k\\) (explicitly T{V]E,k\\) 
TT{V]E,k\\) and N^{V]E,k\\) vNo,{V] E,k\\)). Fur- 
ther, for a visual inspection and comparison of the 
X{V]E^k\\) on the same scale it is convenient to nor- 
malize them by 



x(y-E,k\\) 



X{V;E,kn) 



(6) 



VA^ 

The correlation coefficient, alternatively expressable as 
Ca{V; E) = Aj.^ , is a measure for the 'overlap' of T and 

Na or of T and N^. (i) Consider a constant transmit- 
tance and a constant spectral density, T{V; E,k\\) = t 
and Na{V; E,k\\) = n^, which can be viewed as 'com- 
pletely overlapping'. Then Aj^2 = t^, A^g = n^, and 
At No, = tUa- Consequently, = I which we will con- 
sider as perfect correlation, (ii) Consider a T{V; E^k\\) 
which is nonzero only in a region Qt of the 2BZ: 

T(T/;S,fc||) = T(T/;S,fc||)lo,(fe||), (7) 

with the indicator 



i^t(^ii) 



1 k\\ e Qt 
ku Qt 



(8) 



Likewise, N^iV; E^k\\) is assumed nonzero in a region 
r^Ar« which is disjoint with Qt {^t^^No, = 0; 'zero over- 
lap'). Consequently, Arp2 ^ 0, A]y2 ^ 0, and Atn^, = 0, 
giving Ca = 0. We consider this case as perfectly un- 
correlated. We note in passing that there may be other 
definitions of correlation coefficients, as it is the case for 
the correlation of random variables^. 

Finally, the partial density of states PDOS^ = An^ 
and the total density of states TDOS = A^^ 
(= An^) are computed based on N^^ respectively. 
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III. CONDUCTANCE CALCULATIONS 

Because our theoretical approach to spin-dependent 
tunnehng has been described in detail elsewher^^^^^ 
we restrict ourselves to a brief survey here. The electronic 
structure of a tunnel junction is computed within the lo- 
cal spin-density approximation to density-functional the- 
ory, as is formulated in multiple-scattering theory^. Our 
spin-polarized relativistic layer Korringa-Kohn-Rostoker 
(KKR) methocF™ provides on the one hand the Green 
function of the entire system, from which the spectral 
densities Na{V; E^k\\) can be computed, eq. (|3|. On 
the other hand, the Green function can be usea to cal- 
culate the transmission T(V; /cy), eq. ( [2|), wi thin the 
framework of the Landauer-Biittiker theory^^J^. We can 
also follow the approach introduced by MacLaren and 
Butlei'^ in which the Bloch states in the electrodes and 
the scattering matrix of the interface region are com- 
puted by means of layer-KKR algorithms'^. Here, the 
Green function is not computed explicitly. 

The T and N maps have been computed on identical 
k\\ meshes in the entire 2BZ, with at least 40 000 points. 
In the following applications we restrict ourselves to the 
case of zero bias = and /i/: = /i7^). 



IV. APPLICATIONS 
A. Fe(OOl) 

The properties of the proposed correlation analysis 
are best introduced by an — admittedly — trivial case: the 
spin-resolved Sharvin conductance'^ of Fe(OOl). Because 
the three regions of the junction — £, X, and 1Z — are iden- 
tical, the scattering matrix in eq. ([2| is Sxp = 5xp and 
T{k\\) is an integer. Accordingly, the N map is a pro- 
jection of the Fermi surface onto the (001) plane. With- 
out spin-orbit coupling, spin is a good quantum number; 
hence, we treat the majority and minority channels sep- 
arately. 

Figure [l] displays the normalized /cy -resolved trans- 
mittance and spectral density maps. According to Ta- 
ble [ij the later are ordered based on their {s, p, d} or 
{Ai, A5, A2+A2/}-like orbital contributions. The assign- 
ment of the angular momentum orbitals to the different 
subsets a is motivated by the irreducible representations 
of the symmetry group. This decomposition is given 
in Table [D 

(i) The map of the total spectral density agrees nicely 
with the associated transmittance map in the majority 
channel (top row). In particular, all features are present 
and the correlation is consequently sizable {c^ = 0.802). 

(ii) Other maps which essentially capture all features 
within the majority channel are the A^d and A^As-iike 
maps. A visual comparison of both with the trans- 
mittance could lead to the conclusion that the former 
matches slightly better than the latter. The correlation 



Table I. Decomposition a of the spectral density with respect 
to angular momentum {£) and orbital (m) quantum numbers 
according to Eq. ([3|. Incrementing i provides a classifica- 
tion by means of s, p, and d orbitals. An alternative de- 
composition can be obtained with respect to the irreducible 
representations of the point group C^v (26 . Because such a 
decomposition holds strictly speaking only for the 2BZ center 
(r, k\\ =0), we refer to 'Ai-like' maps etc. 



a 


i 


m 


orbitals 


Ai-like 


0, 1, 2 







As-like 


1, 2 


-1, 1 


Pcc, P^y, dsxz, dsyz 


A2-like 


2 


-2 




A2/-like 


2 


2 


d-Sxy 



analysis, however, shows that this may be a misinter- 
pretation; both coefficients indicate sizeable correlations 
(cd = 0.798 and CAs-iike = 0.828) but that of the As-like 
subset is slightly larger. 

(iii) As a result of the relatively small correlation coeffi- 
cients of s (cs = 0.365) and p (cp = 0.432) we identify the 
d majority states as the dominating conducting channels. 
Further, with the help of Table [l| and the Ca of the A-like 
maps, a hierarchy of conducting d states can be speci- 
fied. Since CAi-iike = 0.484 is quite small, the d^z'^_r2 
states seem to play no significant role. Due to the large 
CAs-iike the dsxz and d^yz orbitals appear to form the 
leading transport channels, followed by the dj;2_y2 and 
dxy states (c(A2+A'2)-iike = 0.656). 

(iv) Similar observations can be made for the minor- 
ity channel (bottom row). Again, the d states represent 
the main conducting channels. But in comparison to the 
majority channel, the role of {dsxz^ ^Syz) and (da^2_^2, 
^3xy) orbitals is interchanged. 

This example shows that the correlation analysis pro- 
vides a powerful analysis tool which on one hand fits 
nicely to the visual interpretation of T and N maps. On 
the other hand, it clearly reveals possible misinterpreta- 
tions, as has become evident in point (ii). 

With 94% (majority) and 66% (minority) the d states 
constitute by far the main contributions to the total den- 
sity of states (TDOS). In this example the hierarchy of 
correlation coefficients often reflects the hierarchy of par- 
tial contributions to the TDOS (see Fig.[T]). This ordering 
becomes incorrect for instance in the case of As-like and 
(A2 + A2/)-like states in the minority channel. Here, the 
As-like states represent with 64% most of the TDOS, but 
exhibit with a CA^-uke of 0.788 a smaller correlation than 
the (A2+A2/)-like contributions (26%, ca^+a^, = 0.847). 

This last finding indicates already that a one-to-one 
mapping of partial density of states (PDOS) hierarchies 
and correlation coefficients is not viable. In general, 
transmittances depend not only on the number of avail- 
able states but rely also on other conditions like e.g. 
the wavefunction matching at interfaces. Hence, as a 
prototype which exhibits more complicated correlations 
between transport and electronic structures properties. 
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Figure 1. Correlation analysis of the spin-dependent Sharvin conductance of Fe(OOl). Spin-resolved transmittance and spectral 
density maps are shown for the entire 2BZ (top row: majority channel; bottom row: minority channel). Each map is normalized 
according to Eq. (|6|. The Nj: maps are decomposed with respect to the representations shown in Table [l] The correlations Ca 
are given beneath each Na map. Numbers in brackets represent each map's partial contribution to the total density of states 
^Nocl • Logarithmic color scales for T and Noc are given at the bottom. 



Fe/MgO/Fe MTJs will be discussed now. 



B. Fe(OOl) /MgO /Fe(OOl) 

In the following an Fe(001)/MgO/Fe(001) MTJ com- 
prising six monolayers MgO is analysed. The magnetic 
directions within both Fe leads are collinearly aligned 
to each other and considered for the case of a parallel 
magnetic configuration. Below we discuss the correlation 
analysis of the majority channel in more detail, because 
the current is dominated by the spin-up carriers^. 

The corresponding transmittance map is displayed in 
Fig.jij As typical for f(k\\) in the majority channel [8 
a Gaussian-like, radially symmetric distribution is found 
around the 2BZ center. 




Figure 2. Normalized majority transmittance map T(h\\) of 
an Fe(001)/6MgO/Fe(001) MTJ. The color scale for f is 
logarithmic. 



Let us suppose briefly that the Julliere^ model is valid 
and that the transport properties can be interpreted ex- 
clusively based on the available PDOSq,. Then the spin- 
polarized conductances are estimated with the product 
of the densities of states at the Fermi energy, strictly 
speaking with the PDOS^ inside the left (£) and right 
lead (7^). In this picture, the conducting channels would 
be, according to the previous discussion of iron, mainly 
determined by d-like or {A5, A2 + A2/}-like Bloch states 
of the Fe electrodes. But a visual inspection and com- 
parison with the majority maps in Fig. [l] reveals at first 
glance no similarity in the structures with neither the 
T(k\\) nor the A^d(^||) maps. 

Consequently, one could ask which states are essential 
for the transport processes if the predominant Fe bulk 
states do not play a decisive role. A closer look at the A^s 
and A^Ai-iike maps in Fig. [l] leads to the identification 
of centrosymmetric blobs like that in Fig. [2] But the 
associated Bloch states exhibit small Cq, and marginal 
contributions to the TDOS in Fe bulk. Do these states 
play a decisive role for the transport within the MTJ? 

To answer this question we consider in a next step the 
layer-wise contributions of the PDOS^ to the respective 
TDOS at the Fermi energy (see Fig. [3|. In particular, 
the ^N^l percentages are shown for one half of the 
symmetric MTJ and are classified again by means of {s, 
p, d} or {Ai, A5, A2 + A2/}-like orbital decompositions. 

The values for the outermost Fe layers in Fig.[3](a) and 
(b) are identical to those for the pure Fe lead in Fig. [l] 
indicating the bulk-like character of those layers far from 
the interfaces. The corresponding d states represent with 
80%-95% the most numerous parts of the TDOS up to the 



5 





100 




80 




60 








40 








20 







(a) 




s 








r 


m 






c 


I ▲ 





















(b) 

Ai-like ♦ 
Ag-like 
(A2 + A2')- like A 



<f <e <e "^""^^^ 




O^Xpr^ ^ ^ ^ ^ 

000 o 00 



Figure 3. (Color online) Layer- and orbital-resolved con- 
tributions of the partial density of states to the to- 
tal density of states within the majority channel of an 
Fe(001)/6MgO/Fe(001) MTJ. The fractions are shown for 
(a) s, p, d, and (b) {Ai, A5, A2 + A2/}-like orbital decompo- 
sitions. 



second Fe layers adjacent to the MgO interfaces. At these 
interfaces the number of d states is drastically reduced. 
Within the MgO film the decrease is continued down to 
about 10% inside the middle region of the tunnel barrier. 

On the other hand, the s and p contributions which 
are apparently vanishingly small within the Fe electrodes, 
obtain substantial weight inside the MgO spacer. In par- 
ticular, the p states exhibit a maximal percentage of 60% 
within the 1. MgO monolayer. Deeper inside the MgO, 
this fraction reduces to about 40%. In contrast, the s 
fractions reach a level of 20% at the Fe/MgO interfaces 
and increase monotonously to roughly about 40%. Thus, 
in the middle of the MgO the number of p and s states are 
comparably large, with a slight advantage of the former. 

Further, a decomposition in {Ai, A5, A2 + A2/}-like 
contributions reveals completely different characteristics 
in Fig. [sJd. Here, the Ai-like contributions, which are of 
minor importance within the Fe leads, experience a mas- 
sive increase at the Fe/MgO interfaces and reach levels 
of about 80% within the tunnel barrier. The remain- 
ing 20% are predominantly occupied with As-like and 
few (A2 + A2/)-like states. This hierarchical order of 
PDOSa within the MgO reflects the well known fact of 
symmetry-selective decay lengths of the evanescent Bloch 
states within the tunnel barriei*^^. 

Due to their sizable presence within the bottle neck of 
the junction, i. e. the MgO barrier, s and Ai-like states 
might indeed characterize the transport. But whether 
this predominance also results in a dominance of the con- 
ducting channels can only be answered by an analysis of 
the transmittance map in Fig. [2] with the respective 
maps. In order to specify exactly which layer-resolved 
spectral densities fits best the structures of the transmit- 
tance map, the maps of the whole Fe/MgO/Fe tunnel 
junction have to be inspected. These maps are shown in 
Fig. [4^. 

One can see immediately that N^, TVg, ^p, and 
^Ai-iike within the MgO layers exhibit a great similarity 



with the structure of T{k\\) in Fig.[2|^ On the other hand, 

the similarities of A/'^, A'As-iike and A/'(A2+A2/)-iike within 
the same layers appear rather small. Due to the great dis- 
similarities of the spectral density maps within the Fe lay- 
ers it is reasonable to expect low correlations there. But, 
based on visual comparisons it is hard to make defini- 
tive statements about similarities of the structures and 
hence which states provide the largest contributions to 
the electronic transport of the whole MTJ. 

At this point, the tool of the correlation analysis pro- 
vides the potential to gain clearer statements. The com- 
puted correlation coefficients which represent a measure 
for the similarity of two maps are summarized in Fig. [4J3 
and [it. 

(i) Considering Figure [IJd, it turns out that s states 
exhibit the most prominent correlations in all layers. For 
each Fe layer the Cg can be quantified with about 0.1- 
0.2 as nearly double as large as those of p and d states. 
Together with the significantly high correlations of 0.7 
(1. MgO layer) up to 0.95 (3. MgO layer) inside the tunnel 
barrier, the previously assumed dominant role of the s 
states can be regarded as proven for the entire MTJ. 

(ii) However, inside the MgO layers the p and d states 
exhibit sizeable increases of their correlation coefficients, 
too. In particular, the Cp are with 0.2 (1. MgO layer), 
0.4 (2. MgO layer) and 0.7 (3. MgO layer) approximately 
twice as high as the c^. In the discussion of the CAi-iike 
coefficients below it will become evident that the increas- 
ing Cp and Cd are mainly related to states exhibiting Pz 
and d322_^2 orbital character. 

(iii) Since the TDOS maps comprise structures of s 
states closley correlated to the transmission map and less 
correlated p and d states, the correlation values of are 
always lower than those for the s states. 

(iv) Along the atomic layers of the MTJ, the char- 
acteristics of the CAi-iike coefficients in Fig. ^ show a 
qualitatively similar dominance as it was found for the 
Cg coefficients in Fig. Within the Fe electrodes the 
correlations of the A 1 -like states are just slightly but no- 
ticeable larger than those of the s states in these layers, 
indicating an additional relevance of Pz and d^z'^-r'^ or- 
bit als in the electronic transport. This finding is sub- 
stantiated by relatively small correlations of both other 
A-like representations inside the tunnel barrier. Since 
the latter exhibit only {p^^, p^, d^xz^ ^Syz} and {da;2_^2, 
dxy} orbital character (see Table [l|, it is reasonable to 
assume that the rising importance of p and d states in- 
side the MgO is identical to an increasing significance of 
Pz and dsz'^-r'^ orbitals. 

The comparison of the results of the correlation analy- 
sis with the discussion of available PDOSq, in Fig.[3]shows 
common features but reveals also significant differences. 
A common outcome of both approaches is the principal 
importance of the tunnel barrier for the electronic trans- 
port of the MTJ. In particular, both discussions end up 
with a conclusion that Ai-like states — i. e. s, p^ and 
dsx'^_r2 orbitals which are preferentially aligned along the 
transport direction — tunnel most effectively and conse- 
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Figure 4. Correlation analysis for the majority channel of an Fe(001)/6MgO/Fe(001) MTJ in parallel magnetic configuration, 
(a) Layer- and orbital-wise ordered spectral density maps (left right) are shown for the entire 2BZ with decompositions 
OL according to Table [l] (top bottom). Correlation coefficients Ca for contributions (a) with a = s, p, d and (b) ot — Ai-like, 
As-like, (A2 + A2/)-like. Color scales for No. are logarithmic. 



quently carry the dominant part of the tunnel current. 

The fact, that Ai-like states, especially those with s 
orbital character, show their decisive role also within the 
Fe leads, represents a qualitative different outcome of 
the correlation analysis. In contrast, these states exhibit 
the lowest PDOSq, contributions within the Fe layers in 
Fig-H 

In principle, the symmetry- selective filtering of the 
MgO tunnel barrier shows up by means of the hierarchi- 
cal order of the {Ai, A5, A2 + A2/}-like PDOSc, fractions 
in Fig. [3j But deeper within the MgO these percentages 



stay rather constant. The layer- wise increase of the effec- 
tive tunneling processes of Ai-like states inside the MgO 
can only be seen by the increasing characteristics of the 
CAi-iike coefficients in Fig. [4]3. 



V. CONCLUSIONS AND REMARKS 

The analysis of transport properties within contacts 
that exhibit planar geometries is often accompanied with 
a visual comparison of large ensembles of A^y -resolved lo- 
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cal spectral density and transmission maps. In this ar- 
ticle we presented an analysis tool which helps to avoid 
laborious inspections and potential misinterpretations by 
providing exact measures. Applying the proposed cor- 
relation analysis the contributions of atoms or orbitals 
can be quantified unambigously. The dominance of the s 
states in the tunnel current of Fe/MgO/Fe MTJs could 
be proven even for the states in the Fe leads, where the 
contribution of states to the TDOS is only marginal. 

Further, the method extends the popular discussion 
of transport proporties at the f point in Fe/MgO/Fe 
MTJs to a more comprehensive anlysis which comprises 
the entire two-dimensional Brillouin zone. In principle, 



the technique is not restricted to planar junctions. Be- 
side more complex layered structures it can be applied to 
other nano structures like point contacts or nano wires. 
For these atomic-sized contacts, we expect that the corre- 
sponding correlation coefficients provide the same phys- 
ical insights as one would obtain from an analysis of the 
conduction eigenchannels [28]. 
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